Dynamics of shallow impact cratering 
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We present data for the time-dependence of wooden spheres penetrating into a loose non-cohesive 
packing of glass beads. The stopping time is a factor of three longer than the time d/v a needed to 
travel the total penetration distance d at the impact speed v a . The acceleration decreases monoton- 
ically throughout the impact. These kinematics are modelled by a position- and velocity-dependent 
stopping force that is constrained to reproduce prior observations for the scaling of the penetration 
depth with the total drop distance. 
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Granular impact is a phenomenon of natural interest. 
One focus of recent work is the size and morphology of 
the crater, and the analogy with planetary cratering Jl|, 
0, Q ■ Another is the dramatic splash produced by the 
collapsing void 0, Still another is the final depth 
of penetration, d, because it probes granular mechanics 
via the depth-averaged stopping force: (F s ) = mgH/d, 
where m is the projectile mass, g = 980 cm/s 2 , h is the 
free-fall distance, and H = h+d is the total drop distance 
BHHHEil For shallow impact by spheres [a. [Tc3| . and 
for deeper impacts by cylinders with various tip shapes 
0, the penetration depth scales as 



d/d = (H/d a ) 1/3 , 



(1) 



where d is the minimum penetration for h = 0. The in- 
set of Fig.^shows Eq. (JTJ agreeing with data over nearly 
three decades in H . While this constrains the stopping 
force, it does not reveal a unique form. For example, 



Eq. (JJJ is equally consistent with F l 



' s cx z 2 and F s cx v 4 ^ 3 , 
where z and v respectively are the instantaneous depth 
and speed of the projectile |fj. Which, if either, of these 
possible stopping forces is correct? How does the nature 
of the stopping force, and the resulting transfer of energy 
from the projectile to the medium, conspire to produce 
the subsequent granular splash and the final crater mor- 
phology? The unintuitive response of granular media to 
external forcing is a topic of wides prea d interest beyond 
the specific example of impact [ll|, [lj, [lj| . 

Recently impact d yna mics have been measured by 
high-speed video [lj, [l5j, ha and by an embedded ac- 
celerometer E3- In Ref. jlj the total upward force is 
found to be HF = —mg+ (mg + kd) . The stopping force, 
in parentheses, is independent of time but has a value 
that depends on the impact speed v ■ Solution of Y,F = 
ma give s a penetration depth of d = y/ mv Q 2 / (2k) . In 
Ref. [lfj, the total force is found to be SF = —mg + k\z\. 
The stopping force is Coulomb friction and increases 
with time. Solution of YiF = ma gives a penetration 
depth of similar form to Eq. ffl: d/d = (H/do) 1 / 2 with 
d = 2mg/k. In Refs. [luLLD], the acceleration decreases 
with time. The various reported force laws thus appear 



contradictory, both in terms of their time dependencies 
and in terms of their predicted penetration depths. Fur- 
thermore, none of the reported force laws is consistent 
with the d ~ H 1 ^ observation of Eq. Q. 

In this paper we measure cratering dynamics in the un- 
explored regime of shallow impact, where the projectile 
never submerges. Our approach is to measure position 
vs time with an optical method, both faster and more 
precise than imaging. Like Refs. [lfil Il7j. we find that 
the acceleration decreases throughout impact. Our the- 
oretical approach is to consider possible instantaneous 
force laws whose depth-averages reproduce the observed 
scaling of Eq. 0J. The best candidate depends on both 
position and speed, and su gges ts that the seemingly dis- 
parate results of Refs. 0,111 EE EH 03 may not be 
contradictory but instead may represent limiting cases of 
a common force law that holds for both shallow and deep 
impacts. 

Our materials and penetration depth measurements 
are identical to those of Ref. E3 • The medium consists 
of glass beads, diameter D g = 0.30 ± 0.05 mm, prepared 
at 59% packing fraction by slowly turning off a fluidizing 
upflow of air. The projectiles are wooden spheres of di- 
ameter Df, = 1.49 inch or Dt, — 1.99 inch, and of density 
Pb = 0.7 g/cc. These are held and dropped from rest via 
a suction mechanism. The drop distance and final pene- 
tration depth are measured with a telescope mounted to 
a height gauge. Data for the wooden spheres of Ref. [Io| . 
and for the dynamics runs reported here, are shown in 
the inset of Fig. [I] to obey Eq. {Q . 

The time-dependent vertical position of the projectile 
is measured optically. A laser (A = 532 nm) is placed 
about 1.5 m from the sample, along with a cylindrical 
lens that fans the beam into a thin sheet. An aperture 
is used to select the central portion of the beam, where 
the intensity is nearly constant, and to set its size to be 
slightly greater than D^- On the other side of the pro- 
jectile, we align a second aperture equal to D^. Behind 
this we place a large planoconvex lens to focus the light 
onto a photodiode. As such, the collected light intensity 
varies nearly linearly with projectile position. In each run 
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FIG. 1: The vertical position z of the bottom of a sphere vs 
time t, for typical conditions as labelled. The impact begins 
at {z = 0,t — 0}, and gravity points in the — z direction. 
Only every 1/100*' 1 datum is shown. The inset shows the 
final penetration depth vs total drop distance, both scaled 
by the minimum penetration depth d c . The smal l sy mbols 
represent data for four wooden spheres from Ref. 10]. The 
large symbols represent new runs, for which in this manuscript 
we report on impact dynamics (o is for the example in the 
main plot). The values of d are 0.63 cm and 0.44 cm for the 
1.99 inch and 1.49 inch diameter spheres, respectively. 



the photocurrent starts at a maximum, decreases as the 
projectile freely falls onto the beam, goes to zero when 
the beam is fully blocked, and then increases as the pro- 
jectile falls further and light passes over its top; impact 
occurs during this last phase. To calibrate we hold the 
ball in several known positions, with light passing above 
and below, and we fit to the particular cubic polynomial 
expected for an aligned Gaussian beam. 

Typical depth vs time data are displayed in Fig. ^ 
The impact occurs at t = 0, and the projectile comes 
to rest in about 0.03 s. The solid curve through the 
t < free-fall data is not a fit, but rather — v t — gt 2 /2 
with w = \/2gh. We digitized the photocurrent using a 
12-bit A/D converter operating at 10 5 points per second. 
Finer digitization and faster capture could be achieved at 
marginally greater expense. The fidelity with which im- 
pact dynamics can be captured by our method exceeds 
high-speed video. Notwithstanding, there exist certain 
limitations. One is that for large enough drop heights 
the grains splash into the laser beam. A lesser limitation 
is that the minimum drop height was such that none of 
the beam was blocked by the projectile prior to its re- 
lease. Overall we achieved over a factor of two variation 
in total drop distance for each sphere. To differentiate 
position vs time data, we fit to a cubic polynomial with 
Gaussian weighting that nearly vanishes at the edges. For 
fitting windows that are too small, the velocity and accel- 
eration results are noisy; for fitting windows that are not 
too large, the depth-averaged acceleration equals gh/d as 
required by energy conservation. This check gives con- 
fidence in both our data and our differentiation proce- 
dures. 
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FIG. 2: Vertical components of position, velocity, and acceler- 
ation vs time, for all runs; red (orange) dots are for D), = 1.99 
inch (1.49 inch). To within measurement errors, the data col- 
lapse when lengths are scaled by d and times are scaled by 
d/v , where d is the final penetration depth and v is the 
impact speed. By construction, the scaled position data all 
decay from to —1 with an initial slope of —1 (gray lines). 
Note that the acceleration is not constant (gray dashed line 
ad/vo 2 = 1/2), but rather decreases with time. The data 
rule out the power-law forces consistent with Eq. Q, but are 
consistent with a modified Poncelet model. 



Data for position-velocity-acceleration vs time, for 
both spheres and all drop heights, are displayed in Fig.|2j 
All lengths are scaled by the final penetration depth d, 
and all times are scaled by the time d/v required to move 
a distance d at the impact speed v . To within measure- 
ment error, the data all collapse according to this scal- 
ing. By construction, the scaled position data must decay 
from to —1 with an initial slope of —1; the scaled veloc- 
ity data must decay from —1 to 0; and the depth-average 
of the scaled acceleration must equal 1/2. We find that 
the spheres all come to rest at about 3v a t/d, roughly 
three times longer than if they moved the same distance 
at constant speed. Since v ~ h 1 / 2 and d ~ if 1 / 3 have 
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similar scaring, the impact duration is nearly constant 
as in Refs. 0,0. We also find that the acceleration 
decreases with time, in accord with Refs. 0, but 
in contrast to Refs. [lj,ll5|- Note that the scaled value 
of gravitational acceleration is —gd/v 2 = —d/(2h); this 
ranges from —0.05 to —0.10 for our runs (shaded gray 
region in Fig. [5J;) and is generally small compared to the 
projectile acceleration. 

Our impact dynamics data can now be compared with 
expectations for various candidate force laws. For ex- 
ample the simple ad-hoc form XF = —mg + k\z\ a \v\P 
best agrees with Eq. Q if the exponents are related 
by P = (4 — 2a) /Z. The agreement becomes exact for 
£F = —mg + k\z\ 2 . Predictions for this special case, 
{a = 2,P = 0}, and also for {a — 0,/3 = 4/3}, are 
displayed with the scaled data in Fig. [2J Evidently, the 
decay of z(t) is too fast for the former and too slow for the 
latter. The actual behavior lies between these extremes. 
A marginally-acceptable fit, also shown, is attained for 
{a=l/2,P=l}. 

Better fits can be achieved if the stopping force equals a 
constant plus a term that grows with speed. If the drag 
is viscous, then the force law is given by the Bingham 
model, £F = — mg + (F Q + b\v\). If the drag is iner- 
tial, then the force law is given by the Poncelet model, 
XF = —mg + (F Q + cv 2 ). The Bingham model has re- 
cently been advocated for granular impact ||, while the 
Poncelet model has long been used for ballistics appli- 
cations 0. For both, position vs speed can be found 
by writing a = vdv/dz, separating variables, and inte- 
grating. All our dynamics data are shown again in phase 
space plots of scaled velocity and acceleration vs depth, 
as well as acceleration vs speed, in Fig. The best one- 
parameter fits to these models give (F c — mg)/(bv ) — 
0.065 and (F — mg) / (cv 2 ) — 0.20, respectively; x 2 is 
smaller for the Poncelet model by a factor of two. These 
fits (not shown) are both acceptable, but not quite as 
nice as the one displayed. Still, neither model predicts 
the observed penetration depth scaling of Eq. (JJJ. The 
Bingham model gives d = (mv /b)[l — aln(l + 1/a)] 
where a — (F —mg)/(bv ), and the Poncelet model gives 
d = [m/(2c)]ln(l + 1/a) where a = (F — mg)/{cv 2 )- 
Furthermore the predicted initial accelerations, and the 
final penetration depths, are unphysical in the limit 
v -> 0. 

We now modify the Poncelet model along the lines sug- 
gested by Tsimring and Volfson jl9| ■ Since the constant 
term represents friction, it should depend on depth ac- 
cording to hydrostatic pressure [la,l23 an d the geometry 
of the granular medium near the projectile: XF = —mg+ 
[F(z) + cv 2 ]. For any F(z) this can be recast as an ordi- 
nary differential equation for kinetic energy vs position 
and solved by use of an integrating factor, exp(— 2cz/m). 
To exactly recover the observed penetration depth scal- 
ing, Eq. QJ, and hence be consistent with the observed 
depth-averaged stopping force, we find that the friction 
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FIG. 3: Velocity and acceleration vs position, and accelera- 
tion vs velocity (inset), for all runs, scaled as in Fig. H red 
(orange) dots are for Dt — 1.99 inch (1.49 inch). Excellent 
fits to the modified Poncelet model are shown as solid purple 
curves. The various-dashed curves represent the same power- 
law forces as in Fig. 



term must be F(z)/mg = [3(z/d ) 2 -l] exp(-2|z|/rf 1 ) + l 
where d\ = m/c. This vanishes for shallow penetration 
and approaches a constant for deep penetration, consis- 
tent with the special limits advocated in Ref. [19(. Alto- 
gether the projectile acceleration in our model is 

a/g = [3(z/d ) 2 - 1] exp(-2|z|/di) + « 2 /Mi)- (2) 

For high speeds and deep impacts the velocity term dom- 
inates. For shallow impacts, the leading behavior is 
a/g w -l+2\z\/d 1 + (3/dl-2/df)z 2 +v 2 /(gd 1 ). For zero 
impact speed, the initial acceleration is —g, as expected. 

To compare our force model with data, it is convenient 
to work in phase space since velocity vs depth can be 
computed directly as 



1 - 



\z/d\ 3 -\z/d\(d /d) 2 
1 - (d Q /d) 2 



,-2\e\/di 



(3) 



For v — 0, the relation is v 2 — 2gd (\z/d \ — 
\z/d \ 3 ) exp(— 2\z\/di). Since the predicted dynamics de- 
pend on the value of d/d , the apparent collapse of data 
in Fig. |2 would then represent clustering about an aver- 
age to within measurement error. Setting d/d = 2.77, 
the average value for all our runs, the best one-parameter 
fit to Eq. (PJ) is for d/d\ = 1.2. This gives an excellent de- 
scription of the data, as shown by the solid purple curves 
in Figs. EE 
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As the ball comes to rest, in all acceptable fits, the 
time becomes about 3d/ v and the upward acceleration 
becomes one to two g. After stopping the acceleration 
abruptly vanishes, as seen directly in measurement by 
accelerometer ^t|. The discontinuity is more apparent 
here in plots of acceleration vs position or speed (Fig. I3J) 
than in plots of acceleration vs time (Fig.[2t). 

In conclusion, we have measured the dynamics of im- 
pact to further constrain the form of the force law respon- 
sible for prior observations of final penetration depth vs 
total drop distance, d ~ 7J 1 / 3 [1,01 El- For two sphere 
sizes, and a factor of two in drop distance, the position vs 
time data appear to collapse when scaled by impact speed 
and final depth. Several possible force laws can now be 
ruled out altogether. Two velocity-dependent force laws 
are consistent with dynamics but not with penetration 
depth data. The best candidate is a modification of the 
Poncelet model along the lines suggested in Ref. [19(. It 
features inertial drag plus a particular depth-dependent 
friction term designed to exactly recover the penetra- 
tion depth scaling of Eq. The predicted kinematics, 
Eqs. agree very well with our new data. The model 

has two important length scales, the minimum penetra- 
tion depth d and an inertial drag length d\ . The former 
is given by d Q — (0.14/ p) 3 / 2 (p b / p g ) 3 / 4 D b where fi is the 
tangent of the repose angle, p g is the grain density, p& 
is the ball density, and D b is the ball diameter ^(J- F° r 
the system examined here, d\ is about twice as large as 
d but we do not yet know the dependence of d\ on sys- 
tem properties. This could be deduced, and the model 
could be further tested, by measurement of impact dy- 
namics over a broad range of conditions. At one extreme, 
the "dry quick sand" examined in Ref. |l5l | for v a = is 
reproduced very well by our model with d\ = 3d a . 
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